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Abstract 

We examine critical properties of the quarter-filled one-dimensional Hubbard model 
with dimerization and with the onsite and nearest-neighbor Coulomb repulsion U 
and V. By utilizing the bosonization method, it is shown that the system exhibits 
an Ising quantum phase transition from the Mott insulating state to the charge- 
ordered insulating state. It is also shown that the dielectric permittivity exhibits a 
strong enhancement as decreasing temperature with power-law dependence at the 
Ising critical point. 
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The Bechgaard salts (TMTSF) 2 X and the Fabre salts (TMTTF) 2 X form a 
family of organic conductors with quasi-one-dimensional band structure, with 
a long history of both theoretical and experimental research. These materi- 
als have a 3/4- filled electron (i.e., quarter- filled hole) band due to a charge 
transfer from two organic molecules to one counterion X. Since the transfer in- 
tegral is slightly dimerized along the chain axis, the system becomes effectively 
half-filling. Optical and photoemission experiments on the normal states have 
shown evidences of a charge gap for both TMTTF salts and TMTSF salts 
[1,2]. This charge gap is believed to be related to the formation of a uniform 
insulating state. Recently, a charge ordered (CO) insulating state have been 
discovered by NMR measurements in the (TMTTF) 2 AsF 6 and (TMTTF) 2 PF 6 
salts [3] and by dielectric permittivity measurements in the (TMTTF) 2 PF@ 
salt [4]. 
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In order to describe theoretically the insulator, one has to take into account 
the commensurability effect of both half-filling and quarter-filling [5]. As a 
simplified model, the quarter-filled one-dimensional (ID) extended Hubbard 
model with dimerization and with the onsite and nearest-neighbor Coulomb 
repulsion has been proposed and it is shown in a mean-field approach that a 
transition from the state with a uniform charge distribution to the CO state 
takes place by increasing the intersite Coulomb repulsion [6]. In order to obtain 
the characteristics of such a transition, the effect of quantum fluctuations must 
be taken into account. For this purpose, the bosonization method is applied 
to this system [7,8]. The model reduces to the double-frequency sine-Gordon 
(DSG) model for the charge degrees of freedom where two kinds of non-linear 
terms appear due to the half-filled and quarter-filled umklapp scattering. The 
DSG model has been analyzed classically [9], and quantum field theoretically 
[10-12,8]. In terms of the renormalization group (RG) method, it is shown 
that there are two-kinds of different fixed points corresponding to the Mott 
insulator and the CO insulator [8]. However, in the perturbative RG scheme 
the both fixed points are strong coupling where the perturbative RG treatment 
breaks down, thus the low energy states behavior and the critical behavior of 
the transition cannot be captured correctly. 

Recently, it has been established that the critical properties of the DSG model 
is described as a quantum Ising transition, i.e., c = 1/2 CFT [10-12]. By ex- 
panding the previous analysis of the DSG model, we examine the critical 
properties of the quarter-filled ID extended Hubbard model with dimeriza- 
tion, and also examine the physical quantities, such as dielectric response and 
optical conductivity, in the present study. Note that in [11,12], the related 
problem of the conductivity in a ID Hubbard model with a staggered poten- 
tial has been considered. We also discuss the experimental results on the huge 
enhancement of the dielectric permittivity obtained in the (TMTTF) 2 X salts 
[4]- 

We consider a ID extended Hubbard model with quarter-filling, where the 
corresponding Hamiltonian is given by 

H = - E [* + (-!)''**] (4 c ,+i,* + H 
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where n^ a = Cj a Cj tU and c^ a denotes the annihilation operator of the electron 
at the j-th site with spin <t(=|,|), and t d (> 0) denotes dimerization along 
the chain. The onsite and nearest- neighbor Coulomb repulsions are denoted 
by U and V, respectively. By utilizing the bosonization method, the phase 
Hamiltonian is given by H = H c + H s with [8] 
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where 0± and <p± are charge and spin phase variables [8] and a is a cutoff of 
the order of the lattice constant. In Eq. (2), the velocity and the Tomonaga- 
Luttinger exponent are given by v p = v F [(l + g 4p / 2ttv f ) 2 — (g p /27iv F ) 2 } 1 ^ 2 , v a = 
v F [(l+g 4a /27rv F ) 2 -(g a /27rv F ) 2 } 1 / 2 , K p = [(2TTV F +g 4p -g p )/(2nv F +g 4p +g p )} 1 / 2 
and K a = [(2irv F +g 4a -g a )/ (2-Kv F +g 4<T + g a )} 112 . The constants g p , g 4p , g a and 
g 4u are given by g p = g 4p — (U + 4V)a and g a = g ia = Ua in the lowest order 
of the interaction, where a is the lattice constant. The magnitude of half-filled 
umklapp scattering, g 1 / 2 , and that of quarter-filled one, gi/ 4 , are given by [8] 



"'''vkr^v-'^ (4) 

and gi/2 = BUa in the lowest order where A = [1 — (td/t) 2 ]/[l + {td/t) 2 } and 
B = [2td/t]/[l+(td/t) 2 }. We neglect the spin-charge coupling term U 3 cos46* + cos20_| 
due to its large scaling dimension [13]. The order parameters of the Ak F charge 
density wave (CDW) state and the Ak F bond ordered wave (BOW) state are 
given by [13] 



4kF - C Dw(x) oc (-l) x / a cos26 + (x), (5) 
4kF - BO w(x) oc (-l)*/ a sin20 + Or), (6) 



where k F = 7r/4a. 



The spin degrees of freedom becomes gapless due to the marginally irrelevance 
of the g s -term. Therefore we focus hereafter on the charge degrees of freedom 
only. In the case of tj = (i.e., gi/ 2 = 0), it has been shown by the RG analysis 
that the system undergoes the Kosterlitz-Thouless transition by increasing the 
nearest-neighbor repulsion [13] where the gapped phase corresponds to the CO 
insulator. In the case of td ^ 0, the Hamiltonian H c has two possible ground 
states with (0+) = 0, tt/2 (4£; f -CDW, dominated by g 1/4 < 0) or (0+) = tt/4 
(4/c^-BOW, dominated by gi/2)- These two ground states are incompatible, so 
we expect a phase transition of the Ising type, as discussed in [12]. 

By following the method of Ref. [12], we analyze the DSG model Eq. (2). By 
using the Jordan- Wigner transformation, the vertex operators are expressed 
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in terms of the order (<7i and a 2 ) and disorder and /i 2 ) parameters of the 
two Ising models as follows [14,12]: 



sin 20+ = (j* a*, (7) 
cos 20+ = ^, (8) 

cos40 + = i(o? + (7f). (9) 

The charge part of the Hamiltonian if c can be expressed in terms of the 
quantum Ashkin- Teller model which consists of two quantum Ising spin chains 
coupled by a four-spin interaction term [12]: 



if c = #qi [0-1,02] + H at [01, 02] + ^[01,02], (10) 

where 07 are the Pauli matrices of Z-th chain (Z = 1,2). In Eq. (10), Hqi is 
given by 



h qi = - E E[ J ^W^(^ + 1 ) + A <W]> (11) 

i=l,2 n 

where n denotes the lattice site on the respective chain. This Hamiltonian 
describes two independent quantum Ising (QI) chains in a transverse magnetic 
field (A). In Eq. (10), Hat is the Ashkin- Teller interaction term given by 



H AT = J'J2[^ 
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where this term describes a four-spin interaction. It is noted that this inter- 
action term takes a self-dual form, thus this term does not make the system 
non-critical. Finally, Hh is a "magnetic field"-type coupling, given by 



H h = -hJ2° z i(n)(T z 2 (n). (13) 

n 

By using the Jordan- Wigner transformation, the parameters in Eq. (2) and 
those in Eq. (10) can be related as v p = 2Ja[l-(2J' /nJ) 2 ] 1 / 2 , K p = (1/8)[{1 + 
(2J'/7rJ)}/{l - (2J'/nJ)}] 1 / 2 , g 1/2 = 2ah and g 1/4 = -4a(A - J). Note that 
the sine-Gordon model, Eq. (2), with K p = 1/8, g 1 / 2 = and g 1 / 4: 7^ 0, 
is equivalent to a model of two degenerate massive Majorana fermions, i.e., 
two decoupled non-critical Ising models [14], and, when the nonlinear term 
vanishes (#1/4 — > 0), two Ising models become critical simultaneously, i.e., 
critical properties are described by c = 1 CFT. 



4 



The origin of the Ising criticality (c = 1/2 CFT) in the case of g 1 / 2 ^ can 
be captured easily by considering the limit of large h [12]. In this case, the 
Ising spins are locked as (<7i,<72) = (1,1) or (—1,-1) where <7j is the eigen- 
value of of. This fact implies that the "relative" degrees of freedom becomes 
locked (i.e., this freedom is always non-critical) while the "total" degrees of 
freedom asymptotically decouples and the system can become critical where 
the critical properties are described by c = 1/2 CFT. In order to analyze 
the physical quantities on and near the criticality, we change the basis into 
"total"- "relative" representation, where the relative Ising mode is described 
by r = o\ o"2, while the total mode by e.g., a — o\. The operator identities are 
given by [12] 



erf, T* = olo\, 



ofoZ, r x = a x . 



(14) 



In terms of these spin variables, Eq. (10) is rewritten as H c = H a + H T + H a 
where 



H a = Y\_-Jo z (ji)a z (n+l) + J'a x (n)], (15) 

n 

H T = J2[j'r z (n)T z (n + l)- Ar x {n) -hr z (n)], (16) 

n 

H trr = -^2[j<r z (n) a z (n + l)r z {n) r z {n+ 1) 

n 

+ A a x (n)r x (n)]. (17) 
The Ising transition corresponds the situation when the a-degrees of freedom 



go massless at A y h(2J + J') [12]. With no regard to the Ising transition in 
the a-chain, the r-chain remains frozen in a configuration with both (r z ) and 
(t x ) nonzero. Here we introduce the spin variables ji and v which are dual to 
a and r, respectively, by the Kramers- Wannier transformation: fi z (n+ 1/2) = 
EI™ =i cr x (m), fi x (n + 1/2) = a z (n)a z (n + 1), v z {n + 1/2) = \\ n m=1 r x (m) and 
v x (n + 1/2) = T z (n)r z (n + 1). The spin variables fi and v can be related 
with Hi and /i 2 by ji z = (J> x = f-t x , v z = ^\ an( i vX — ^\^ x - The 

expectation value of /i 2 becomes finite ((//) ^ 0) in the case that a-chain 
is in the disordered phase ((<j z ) = 0), while becomes zero ((/i 2 ) = 0) in the 
ordered phase ((a 2 ) ^ 0). The order parameter of 4/cp-CDW is given in terms 
of the spin variable by 04fc F -cDW ^ cos2# + = [i\ /if = fi z , and then one finds 
((9 4fcp _ CDW ) oc (/i 2 ) 7^ in the disorder phase only. Therefore the transition 
from the Mott insulator to the CO insulator corresponds to this order- disorder 
Ising quantum phase transition of the a Ising chain, where the disordered 
phase of the Ising system corresponds to the CO insulating state. Note that 
the quantity {O^-bovj) takes a finite value and shows no anomaly at the 
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transition point, since the order parameter of 4A; F -BOW is given in terms of 
the noncritical Ising field: 4fcF - B ow ^ sin29 + = a\ = r z . 

Now we examine the dielectric permittivity and the optical conductivity at 
finite temperature. In the energy scale higher than the gap in the r-chain, the 
perturbative analysis in terms of renormalization group method works well [8]. 
The complex conductivity ct(cj) = a\{uj) + i(T 2 (w) (ci = Recr and a 2 = Imcr) 
can be calculated from cr(u) = i\x(w) — x(0)]/^5 where x is the retarded 
current-current correlation function: x(^) = —i J dx J °° dt(\j(x, t),j(0, 0)]) e tujt ~ 
In terms of the complex conductivity, the dielectric permittivity is given by 
e{uj) = ei(u) + ie 2 (c<j) = 1 + (Am/u) a(u), where e\ = Ree and e 2 = Ime. 
By utilizing the RG method, it can be shown that the resistivity p(T) and 
the dielectric permittivity £i(T) behaves as p(T) = l/a\(T) oc T 16Kp ~ 3 and 
ei(T) oc T 2 ~ 16Kp in the case that the quarter-filled umklapp scattering is dom- 
inant [8]. This analysis, however, only holds in the high temperature region. In 
the energy scale lower than the gap in the r-chain, we analyze these physical 
quantities in the quantum Ising model scheme. By using the spin variable p, 
the current operator are given by [12] 



j(x) = --d t 6 + ^-Cd tf i(x), (18) 

71 

where fj,(x) is a continuum variable of p z {n) and and C(oc (r 2 )) is a non- 
universal numerical constant. Here we have used the fact that the r-field is 
always non-critical and has nonzero expectation value. We introduce the re- 
sponse function of the Ising disorder parameter: 



oo 

D {R \k,u) = -i jdx J dte~ lkx+lult (\p.(x,t), /i(0,0)]). (19) 

o 

By using Eq. (18), the current-current correlation function is given by xi 00 ) = 
C 2 u 2 D( R \0,cu). Then the real and imaginary parts of the conductivity o\ and 
cr 2 are given by u x {u) = -C 2 u ImD^(0, u) and a 2 (co) = C 2 ooReD^(0,uj). 
On the Ising criticality, it is known that the response function is given by 
-Dcrit(r-) oc 1/r 1 / 4 where r = (x 2 + r 2 ) 1 / 2 with the velocity being set to unity. 
This power-law behavior is still valid away from the criticality in the short 
range region r 1/A CT , where A CT is the gap in the a-part of the Ising chain. 
By using this asymptotic behavior, the retarded response functions is given 
by [14] 



^2(0, -) = ~0^K B2 ( A 2A) , (20) 

where S = u/AttT, the conformal dimension is A(= A) = 1/16 and B(x,y) is 
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Fig. 1. The temperature dependence of the dielectric permittivity e\(uj,T) for 
oj/W = 0, 0.0005 and 0.001. The inset shows the frequency dependence of ei(u,T) 
with fixed temperature T/W = 0.001, 0.005 and 0.01. 

the beta function given by B(x,y) = T(x)T(y)/T(x + y). The temperature de- 
pendence of the dielectric permittivity €i(u,T) with zero and finite frequency 
is shown in Fig. 1, where W denotes the cutoff of the order of the bandwidth 
and for simplicity we have set C — 1. The frequency dependence of ei(iv, T) is 
also shown in the inset. The dielectric permittivity shows strong temperature 
dependence. In high temperature limit, the dielectric permittivity ex becomes 
unity independent of uj, but as decreasing temperature the huge enhancement 
can be obtained and e\ takes a maximum at T uj. At low temperature 
(T<u), €i decreases and changes its sign, and finally saturates to a value 
that depends on uj. The optical conductivity also increases as decreasing tem- 
perature, and exhibits a maximum at T uj. The limiting behaviors of e\ and 
(7i are given by e 1 {u,T) = 1 + const. /T 7 / 4 and <Ti(w,T) oc uj 2 /T u / 4 for high 
temperature (T 3> uj), and ei(uj,T) oc — I/uj 1 ^ and ax(uj,T) oc l/c<// 4 for low 
temperature (T <^ a;). 

Finally we discuss the experimental observation of the (TMTTF) 2 X salts. 
Recent analysis of the dielectric response exhibits a strong temperature de- 
pendence and a huge enhancement of the dielectric permittivity as decreasing 
temperature [4]. In our analysis of the temperature dependence of the dielec- 
tric permittivity, we have restricted ourselves to the case on the Ising critical 
point. However, the results would be valid in the vicinity of the transition 
point in the energy region above A CT . Therefore the huge enhancement could 
be interpreted by considering that the system is in the vicinity of the transition 
point from the Mott insulator to the CO insulator. 
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In summary, we examined the properties of the dimerized extended Hub- 
bard model at quarter-filling in one dimension. It was shown that the phase 
transition from Mott insulating state to the CO insulating state in this sys- 
tem corresponds to the Ising quantum phase transition described by c = 1/2 
CFT. We have also show that the dielectric permittivity exhibits a strong en- 
hancement as decreasing temperature with power-law dependence at the Ising 
critical point. 

One of the authors (M.T.) thanks A. Furusaki, T. Hikihara, H. Seo, Y. Suzu- 
mura and H. Yoshioka for helpful discussions. 
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